Interaction- Assisted Gap Opening in Moire Graphene Superlattices 
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We show that coupling of Dirac particles in graphene to moire superlattice potential in G/hBN 
heterostructures is strongly affected by interactions. Interacting Dirac particles respond strongly 
to a pseudospin-dependent ('colored') potential arising due to the A-B sublattice modulation and 
relatively more weakly to a scalar potential which is sublattice-blind. As a result, the effective 
potential seen by electrons in graphene becomes predominantly 'colored,' leading to pseudospin- 
dependent Bragg scattering. Many-body effects generate giant enhancement of the gap at the Dirac 
point, which can overtake the enhancement of the minigaps opened by Bragg scattering at the 
superlattice Brillouin zone edges. We predict a power-law dependence for the gap vs. moire pattern 
wavelength, which can be used to directly test the effects of interactions in experiment. 



Electronic states in heterostructures comprised of 
graphene and hexagonal boron nitride (hBN) have inter- 
esting and potentially useful properties p~r[3]. Both mate- 
rials have hexagonal lattice structures with near-pefectly 
matched periods. Placing graphene on hBN, which is a 
wide-bandgap electric insulator, yields structures with 
low disorder, high carrier mobility, and improved de- 
vice performance. A small period mismatch (w 1.8%) 
of graphene and hBN leads to moire superlattices with 
spatial periods reaching Ao ~ 14 nm for perfectly aligned 
crystals. The electronic states in graphene have compa- 
rable deBroglie wavelength and are fully exposed to the 
hBN superlattice potental. A number of interesting phe- 
nomena were predicted in this system: Bloch minibands 
and changes in Fermi surface topology 0HS], replica Dirac 
points 71, and most recently, Hofstadter-type commen- 
surability effects in magnetic field [51 [3], all successfully 
described by a noninteracting Dirac-Bloch model. 

As we argue below, the noninteracting single-particle 
Dirac-Bloch picture, despite all the appeal it may have, 
is not fully adequate. We will see that this picture misses 
some important physics occurring near the Dirac point, 
which results from electron-electron interactions. Specif- 
ically, as we will discuss in detail below, interacting Dirac 
particles respond in a very different way to a scalar exter- 
nal potential which is sublattice-blind and a pseudospin- 
dependent ('colored') potential arising due to the A-B 
sublattice modulation. In the first case, interactions 
generate polarization that screens the potential. In the 
second case, interactions generate sublattice correlations 
that amplify the potential. 

The interaction-induced amplification of colored po- 
tential leads to pseudospin-dependent Bragg scattering 
which has dramatic consequences for electronic states. 
First, a gap at the Dirac point opens up. Second, because 
Dirac mass term is a relevant perturbation in the renor- 
malization group (RG) sense, this gap undergoes giant 
interaction- induced enhancement. Counter- intuitively, 
this enhancement can be much greater than the enhance- 
ment of the side gaps opened by Bragg scattering at prin- 
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FIG. 1: Interaction-induced enhancement of couplings to the 
moire superlattice potential. RG flow is shown for the gap 
at the Dirac point, Ao, and the gap opened by Bragg scat- 
tering at the principal superlattice harmonics, Ai (red and 
blue lines). Renormalization proceeds in two stages described 
in the text: a < A < Ao (stage 1) and A > Ao (stage 2), 
where a = 1.42 A is carbon spacing, Ao is the moire super- 
lattice wavelength. Both gaps are enhanced from the nonin- 
teracting values Ag ', A'°\ however the gap Ao undergoes a 
much larger enhancement and, despite a small initial value, 
eventually overtakes Ai. Both Ao and Ai grow faster than 
the carrier velocity v (green line). Inset shows the density of 
states near the Dirac point with the gap Ai opened at the 
superlattice Bragg energy eo = |?w|b|, |b| = 27r/A . Param- 
eters used: superlattice period Ao = 14 nm, scaling exponents 
P = 0.6, f} v = 0.3, initial values A 0) « 0.02 A™ (see text). 

cipal superlattice harmonics, see FigJT] 

Due to spatial oscillations of the coupling to moire su- 
perlattice, the RG flow for these couplings proceeds in 
two separate stages: stage 1 describing renormalization 
at lengthscales up to the moire wavelength, A < Ao, and 
stage 2 describing lengthscales A > Ao- The lengthscale 
at which RG terminates is controlled by the screening 
length which is set by the distance to the gates when 
those are metallic, as in Ref.[T], or by the screening 
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length in the gate if the latter is realized by a proxi- 
mal graphene layer, as in Ref.[13]. We find that stage 1 
acts as a 'booster,' generating an enhancement of energy 
gaps which is much greater than that of carrier velocity. 
At stage 2 the growth of the side gap Ai stalls whereas 
the gap A continues to grow. As illustrated in Figjl] 
the gap Ao eventually can overtake Ai even if the latter 
starts with a larger microscopic bare value. 

Interaction-assisted enhancement of spectral gaps gen- 
erated by microscopic interactions was proposed for 
a variety of different systems, e.g. for chiral metal- 
lic nanotubes[10] and for Quantum Spin-Hall effect 
in graphene where a gap is opened by spin-orbit 
coupling [llj. In Ref.|12j. renormalization of a gap opened 
by sublattice modulation was analyzed using 1/N expan- 
sion, and a scaling exponent f3 = 16/ir 2 N was found. 
Our results for Ao enhancement are consistent with 
Ref.pJ] for A > Ao, however we find a much stronger 
enhancement at short lengthscalcs, A < Ao, where rapid 
RG growth occurs. This is so because the gap A de- 
pends on a combination of microscopic parameters which 
renormalize independently, see Eq.Q, leading to a much 
steeper growth on stage 1 than on stage 2. 

Besides establishing the hierarchy of energy scales for 
different gaps, our approach may help to shed new light 
on the recently observed insulating state, which arises 
despite the extreme cleanness of G/hBN systems [T3l IT4]. 
This state is observed at charge neutrality (CN) when 
the long-range disorder due to charge puddles is screened 
by gates. These observations, which are not under- 
stood, are particularly suprising because previous mea- 
surements have reported record-high mobilities in this 
system away from CN. Several scenarios are at the ta- 
ble: Ref.[T2] attributes the resistive behavior at CN to 
weak localization due to remnant short-range disorder, 
a picture supported by a 'disorder by order' model [T5] , 
Ref.[14] observes a highly resistive state, with strong in- 
sulating behavior turning on abruptly as a function of 
doping through CN and no appreciable sign of weak lo- 
calization away from CN. This suggests an energy gap 
opening as a source of the insulating behavior. Indeed, 
as pointed out in Ref.[Hj, a gap at the Dirac point may 
be induced by sublattice modulation in hBN substrate. 
Yet without a sizeable boost by interaction effects the gap 
would be too small to explain the observations. Ref. |14j . 
while leaning towards transport gap explanation, did not 
reach a definite conclusion about the origin of the ob- 
served behavior. Presently it is not known which of the 
two scenarios, a bandgap or a transport gap, is responsi- 
ble for the new resistive state. 

Turning to the analysis, since moire superlattice period 
is much larger than the lattice constant, Ao 3> a, states in 
valleys K and K' are effectively decoupled. Hence we can 
describe each valley by an effective continuum Hamilto- 
nian having a nontrivial pseudospin structure reflecting 
A-B sublattice modulation, and an oscillatory position 



dependence reflecting the superlattice periodicity: 
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/(x)=2 ]T cos(b s x + s ), (2) 

s=l,2,3 

P± = Px^Py Here b s (s — 1, 2, 3) are three Bravais vec- 
tors of the triangular superlattice oriented at 60° angles 
relative to each other. Introducing pseudospin Pauli ma- 
trices 0i,2,3 and adding long-range interactions we write 

H = I rf 2 %y^ V4( x ) i vcr P + m 3(x)t73 + ra (x)] ipj(x) 
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where m 3 = §(« u - u 22 )/(x), m = |(u u + u 22 )/(x). 
Here N = 4 is the number of spin/valley flavors, and 
ri(x) = J2i=i n V'l( x )V'i( x ) is particle density. The off- 
diagonal terms u\2 and m 2 i, which describe a gauge-field 
coupling generated by strain, can be incorporated in the 
vap term. However, ab initio studies 17 indicate that it 
is a small contribution compared to un and U2i- 

The amplitude of sublattice modulation can be inferred 
from ab initio calculations JH] predicting 6(iiu — u 22 ) « 
53meV for the equal-period case, b = 0. This gives 
amplitudes 7773 of individual Bragg harmonics which are 
more than 20 times smaller than the kinetic energy at 
the Bragg vector, e = |^ u |b| ~ 150 meV estimated for 
largest superlattice period Ao = 14 nm [TH]. Hence we 
can employ perturbation theory in a small ratio 7713/eo. 

The RG analysis of this Hamiltonian predicts that the 
7773 harmonics grow under RG, whereas the mo harmonics 
do not grow (see below). Therefore, even if the micro- 

copic values m 3 °^ and mj ' are comparable, the bandgaps 
will be dominated by the 777.3 harmonics. We can there- 
fore estimate the bandgap opening at the edge of the 
superlattice Brillouin zone as Ai = 27773. 

In addition, importantly, the interaction with super- 
lattice leads to gap opening at the Dirac point. While 



a single 777,3 harmonic is sign-changing, 



zbx\ 



= 0, and 



thus it cannot give rise to such a gap by itself, a combina- 
tion of three different harmonics can open up a gap at the 
Dirac point. Choosing triplets of harmonics with the sum 
of their wave- vectors adding up to zero, b^ + bj +bfc = 0, 
third-order perturbation theory in 7713 yields a constant 
sublattice-asymmetric term [26 
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where A = 127773/(7; |b|) 2 (the corresponding self-energy 
is shown in Fig[2^,). Contributions similar to Hq can 
also arise when two out of three 777,30-3 terms are replaced 
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FIG. 2: Self-energy describing gap opening at the Dirac point 
(a) and log-divergent diagrams contributing to gap renormal- 
ization at one loop (b,c,d). The contributions (b) and (c), 
accounting for vertex renormalization and velocity renormal- 
ization, arise from integration over |b| < |p| < po (stage 1), 
giving Eq.( |10[ ). The contribution (d) arises from |p| < b 
(stage 2), giving Eq.flTl. 



by too harmonics. However, since these harmonics do 
not grow under RG (see below), these contributions are 
small. Since m-3 <C en, the predicted numerical value 
is small, Ao <C Aj (similar observation was made in 
Ref.[5D]). However, as we find shortly, the gap A un- 
dergoes a giant boost due to interaction effects, growing 
faster than Ai. As a result, the physical values Ao and 
Ai become comparable. 

The effects of interactions can be treated by an RG 
approach, constructed by extending the framework de- 
veloped in Refs. [211 - 55] to account for renormalization of 
the terms too and TO3 in Eq.(|3). There are two distinct 
flavors of RG, a weak-coupling approach using e 2 /hv <C 1 
as an expansion parameter, and a large- N approach, both 
leading to qualitatively similar results [23]. We present 
detailed results for the large- AT approach. 

We first analyze renormalization of the mo and to 3 
terms in Eq.Q treating them as spatially uniform, ig- 
noring the x-dependence. The approach is valid over the 
range of lengthscales a < A < Aq. Larger lengthscales, 
A ^ Ao, will be analyzed below. Renormalization is found 
by dressing the mo and to 3 vertices in Feynman diagrams 
with vertex corrections and analyzing the corresponding 
log-divergent contributions. The vertex correction for the 
Too term is cancelled by a corresponding contribution to 
the Greens function residue, owing to the Ward identity 
that follows from gauge invariance. This can be seen 
more explicitly by analyzing the self-energy 
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with G(e,p) = l/(ie — vap — m^a 3 — mo) the Greens 
function, and V representing the RPA interaction, 



y(e ' p) =i-AV (p)n (e ,p)' ^ o(p) 



«P 



(6) 



Renormalization of the mo coupling, at linear order in 



too, is described by the quantity 9, e S + 3„ ifl S which van- 
ishes due to the form of G(e, p). 

We note that the cancellation of log-divergent contri- 
butions due to Ward identity does not apply to reducible 
diagrams. The latter generate a vertex correction given 
by the RPA dielectric constant e = 1 — NVo(p)H(e, p), 
which describes intrinsic screening of mo vertex by inter- 
band and intra-band polarization. 

Since TO3 vertex is distinct from mo vertex, we expect 
that log-divergent contributions do not cancel. As a re- 
sult, electrons become 'colored,' i.e. their coupling to 
moire superlattice potential is dominated by pseudospin- 
dependent interactions. Renormalization of 7713 was an- 
alyzed in Ref.[12] giving scaling exponent (3 — 16/w 2 N 
which is two times larger than the value fi v = 8/ir 2 N 
found for velocity renormalization in Ref . 23 . This leads 
to RG flow equations for the TO3 coupling and velocity: 
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In & is the RG time parameter, with po 



where £ 

27r/a the UV cutoff. Interestingly, the relation f3 = 2/3 v 
also holds in the weak-coupling approach. The RG flow, 
Eq.|7|, predicts power-law enhancement to the TO3 har- 
monic and velocity for a < A < Aq: 
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Similarly, renormalization of Ao oc m^/v 2 is obtained by 
adding the contributions shown in Fig[2jb,c): 
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where we integrate over |p| < p' < p a . Taking N to 
be large, we find the value 32/ir 2 N which is two times 
larger than the exponent /? found in Ref.[T2] and four 
times larger than the exponent /3 V . As a result, Ao grows 
under RG much faster than Ai: 



12to| _ . ,30-2p v (0) 
A ° - ^2|b|2 - ( A / a ) A ■ 



(10) 



A ten-fold increase in Ai transaltes into a hundred-fold 
increase in A , see Fig{T] Thus, despite a handicap due 
to a small initial value, Eq.Q, the physical values for Ao 
and Ai are in the same ballpark. 

For the sake of generality, and acknowledging an ap- 
proximate character of the scaling dimensions obtained 
from one-loop RG, we shall leave j3 and f3 v unspecified 
in the analytic expressions. An attempt to experimen- 
tally determine scaling exponents was made recently in 
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Ref.pS], where a systematic change of the period of quan- 
tum oscillations with carrier density was interpreted in 
terms of Fermi velocity renormalization, giving a value 
p v = 0.5-0.55. This value is considerably larger than the 
one-loop RG result, /3 V = -f^ s» 0.2. This discrepancy 
is not yet understood. For the plots in Fig[TJ we used 
conservative values /3 V — 0.3, j3 — 2f3 v = 0.6. 

At lengthscales larger than Ao renormalization is sup- 
pressed by spatial oscillations in the m 3 term. In this 
regime m 3 stops flowing and the gap Ai stays con- 
stant, as shown by the horizontal line in Fig [I] (stage 
2). In contrast, because the effective Hamiltonian H' is 
x-independent, modes with |p| < |^ continue to provide 
an enhancement to the gap at the Dirac point, giving 



A (A> A ) = (A/Ac/A (Ao) 



(11) 



The RG flow of A at stage 2 is described by scaling 
exponent j3 which is smaller than the value 3/3— 2f3 v found 
for stage 1. Nevertheless, the growth of Ao at stage 2 is 
still faster than that of velocity. As illustrated in Fig.l, 
A can grow by several orders of magnitude, reaching Ai 
for large enough A. For the parameters chosen in FigjT] 
A (red line) overtakes Ai (blue line), eventually making 
the gap at the Dirac point the largest gap in the system. 

Renormalization of all the parameters in the system, 
including velocity v and the gap A , terminates either 
at an effective screening length set by proximal gate or 
at a lengthscale generated self-consistently through a gap 
opening at the Dirac point. In the first case, the effec- 
tive screening length is determined by the distance to the 
gate or by the screening length for the gate, whichever 
is larger. In the second case, realized for systems with 
remote gates or when the screening length is very large, 
RG terminates at a self-consistently defined lengthscale 
controlled by the gap opened at the Dirac point, 



A* = 271-/2*, hq* = A (A*)/«(A*), 



(12) 



where the A dependence is obtained from RG flow. A 
similar approach was employed in Rcf.[10] to estimate 
interaction-enhanced gaps opened at the Dirac point in 
chiral metallic carbon nanotubes. Under realistic condi- 
tions, as discussed above, we expect A <C eo = ^w|b|. 
Hence the selfconsistent lengthscale satisfies A* 3> Ao- In 
this case, plugging the dependence from Eqs.([8]),( 11 1, we 
find (A*/A ) 1+/3 -^ = to(A )/[A Ao(Ao)] ) giving 



A (A*) = 



hv(X ) 
A A (A ) / 

/3(l+3/3-3|3„) 
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7 = 3/3-2/3,- p' 1 ^' 1 . The predicted power-law 
dependence Aq vs. moire wavelength Aq, Eq.(13l, can be 



used to directly probe the effects of interactions. Setting 
/3 = 2/3„ (see discussion above) , the expression for scaling 
exponent simplifies as: 



7 = 4/3, 



2&(l + 3/3„) /3(2-/3) 



1+Pv 
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for which our one-loop RG gives /3 = 7 w 0.27. 

The predicted scaling Ao oc (Ao/a) 7 can be tested by 
comparing the gap values measured in G/hBN systems 
with different misalignment angles, similar to how it was 
done in Ref.[Jj5] for replica Dirac peaks. It is beneficial 
to be in the regime of small misalignment angles, since in 
this case the scatter in microscopic parameters determin- 
ing bare A °^ values is minimal, and at the same time the 
ratio Ao/a is largest, which maximizes renormalization 
effects. This is the case for structures with near-perfect 
crystal axes alignment, such as those used in Ref. [9 . 

In summary, we predict that interactions have dra- 
matic effect on the electronic properties of G/hBN sys- 
tems. The bandgaps opened up by coupling to sub- 
strate are strongly enhanced by interactions. The pre- 
dicted physical gap values can be many times greater 
than the bare gaps opened due to hBN sublattice mod- 
ulation. The effect is most dramatic for the gap at the 
Dirac point which undergoes a giant enhancement that 
boosts it to a value comparable to or exceeding that for 
principal Bragg harmonics. These effects can be directly 
probed by measuring gap values for systems with differ- 
ent superlattice wavelengths, predicted to obey a power- 
law scaling A oc Aq ■ It can also be probed by doping the 
system, since the gap enhancement due to screened in- 
teractions is considerably smaller than that in undoped 
system where interactions are unscreened. This leads 
to the gap collapsing upon doping, resulting in an en- 
hancement in compressibility when the system is doped 
slightly above the gap. The strong effect on compress- 
ibility makes the interaction physics easily detectable in 
experiment. 
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